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ABSTRACT 

Context. The prominent blazar 3C 279 is known for its large-amplitude variability throughout the electromagnetic spectrum and its 
often y-ray-dominated spectral energy distribution. However, the characterization of its broadband spectral variability still lacks a 
consistent picture, and the origin of its high-energy emission is still unclear. 

Aims. We intend to characterize the spectral energy distribution and spectral variability of 3C 279 in its optical high state. 
Methods. Prompted by an optical high state of 3C 279, we organized an extensive multiwavelength campaign with coverage from 
radio to hard X-ray energies. The core components of the campaign were INTEGRAL and Chandra ToO observations in January 
2006, augmented by X-ray data from Swift and RXTE as well as radio through optical coverage. 

Results. The blazar was observed at a moderately high optical state. A well-covered multifrequency spectrum from radio to hard X- 
ray energies could be derived. During the flare, the radio spectrum was inverted, with a prominent spectral peak near 100 GHz, which 
propagated in time toward lower frequencies. The spectral energy distribution (SED) shows the typical two-bump shape, the signature 
of non-thermal emission from a relativistic jet. As a result of the long exposure times of INTEGRAL and Chandra, the high-energy 
spectrum (0.3 - 100 keV) was precisely measured, showing — for the first time — a possible downward curvature. A comparison 
of this SED from 2006 to the one observed in 2003, also centered on an INTEGRAL observation, but during an optical low-state, 
revealed the surprising fact that — despite a significant change of the high-frequency synchrotron emission (near-IR/optical/UV) — 
the low-energy end of the high-energy component (X-ray energies) remained virtually unchanged compared to 2003. 
Conclusions. Our results prove that the two emission components do not vary simultaneously. This provides strong constraints on 
the modeling of the overall emission of 3C 279. When interpreted with a steady-state leptonic model, the variability among the SEDs 
displaying almost identical X-ray spectra at low flux levels, but drastically different IR/optical/UV fluxes, can be reproduced by a 
change solely of the low-energy cutoff of the relativistic electron spectrum. In an internal shock model for blazar emission, such a 
change could be achieved through a varying relative Lorentz factor of colliding shells producing internal shocks in the jet, and/or 
the efficiency of generating turbulent magnetic fields (e.g., through the Weibel instability) needed for efficient energy transfer from 
protons to electrons behind the shock. 
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1. Introduction at y-ray energies, the majority by the Energetic Gamma-Ray 

Experiment Telescope (EGRET) at energies above ~100 MeV 
Blazars are the most variable and violent type of active galac- (H art man et al. (1999)1 ). It is generally accepted that the radio to 
tic nuclei (AGN). They show rapid variability in all wave- opt i ca i /UV continuum of blazars is synchrotron radiation gener- 
length bands, polarized emission in the radio and optical ate d by relativistic electrons in a magnetized jet. For the high- 
bands, and often superluminal motion. The blazar subclass energy continuum from X - to r -ray energies, two fundamen- 
of AGN combines flat-spectrum radio quasars (FSRQs) and ^ different approaches are being cons i der ed (for a recent re- 
BL Lac objects. Accordi ng to the unified model of AGN (e.g., yiew see> e .g.,| Bottcher (2007)1 >. If protons are not accelerated to 



Urry&Padovani (1995)| , blazars are sources which expell jets ultra . high energies, exceeding the threshold for photo-pion pro- 



close to our line of sight. The discovery by the experiments duction Qn the co . spatially produce d synchrotron radiation field, 

aboard the Compton Gamma-Ray Observatory (CGRO) that high . energy emission might be dominated by Comptonization 

blazars can radiate a large — sometimes even the major — frac- of soft photonfj by the non . therma i jet electrons (leptonic mod- 

tion of their luminosity at y-ray energies marked a milestone els; | M araschi et al. (1992)1 |Dermer & Schlickeiser (1993)| 

in our knowledge of blazars. During the CGRO mission about | Si koraet al. (1994)) . Otherwise, the high-energy radiative out- 

90 blazars were detected by the different CGRO experiments put be dominated by pro ton-synchrotron radiation as well 

; as synchrotron and Compton emission from the decay products 

Send ofpnnt requests to: W. Collmar of pions produce d in photo-pion production processes (hadronic 
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models; e.g., Mannheim & Biermann (1992), Aharonian (2000) 
IMiicke et al. (2003)) . 

The optically violently variable (OVV) quasar 3C 279, 
located at a redshift of 0.536, is one of the most promi- 
nent representatives of this source class. The blazar shows 
rapid variability in all wavelength bands, polarized emission 
in the radio and optical, superluminal motion, and a com- 
pact radio core with a flat radio spectrum. The quasar 3C 279 
was found to be a persistent — though highly variable — 
emitter of y-rays above 100 MeV by EGRET. Subsequently 
the source was frequently observed simultaneously at radio 
and optical bands, and was subject of many multiwavelength 
campaigns during the CGRO era (e.g., Marasc hi et al. (1994)| 



|Hartman et al. (1996)[|Wehrle et al. (1998)| >. A complete compi- 
lation of all simultaneous spectral energy distributions (SEDs) 
of 3C 279 collected during the lifetime of CGRO, including 
modeling of those SEDs with a leptonic jet model, is presented 
in |Hartman et al. (2001a)| The conclusion was that the high- 
energy emission most likely consists — especially during flar- 
ing epochs — of several emission components that vary in- 
dependently. The synchrotron self-Compton (SSC) mechanism 
( [Maraschi et al. (1992)1 |Bloom & Marscher (1996)| l dominates 
the X- and soft y-ray part while at higher y-ray energies an exter- 
nal Compton component (e.g., |Dermer & Schlickeiser (1993)1 
|Sikora et al. (1994)) becomes apparent. 

Despite the enormous observational and theoretical multifre- 
quency efforts, the nature and origin of the high-energy (X-rays 
to y-rays) emission of 3C 279 still remains unclear. Also, the 
connection of the high-energy radiation to the low-energy ra- 
diation is not understood. Cross-correlations between different 
energy bands, particularly between optical, X-ray, and y-rays 
were investigated by Hartman et al. (2001b)] but no consistent 
trends could be found. This poses a serious challenge to the 
quite often adopted single-zone leptonic jet models, in which the 
same electron population produces both the low-frequency emis- 
sion through synchrotron radiation and the high-energy emission 
through Compton scattering. 

In order to improve on this general situation, 
|Collmar et al. (2004)| continued to observe 3C 279 in a 
multifrequency campaign organized around INTEGRAL 
(International Gamma-Ray Astrophysics Laboratory) high- 
energy observations. INTEGRAL was launched in October 
2002. In a campaign in June 2003, 3C 279 was detected 
significantly at hard X-ray energies for the first time by the 
INTEGRAL/IBIS experiment ( [Collmar et al. (2004)) . Those 
high-energy observations were supplemented in X-rays by a 
short Chandra pointing, and by ground based monitoring from 
radio to optical bands. In June 2003, the blazar was found in 
the faintest optical brightness (optical R-band: ~17 mag) of the 
last 10 years, roughly 5 mag fainter than the maximum, and 
about 2.5 to 3 mag fainter than average. That campaign allowed 
|Collmar it al. (2004)] to compile a contemporaneous spectral 
energy distribution (SED) of an exceptional optical low-state. 

In order to measure emission changes as function of op- 
tical brightness, we proposed for an intensive multifrequency 
campaign during an optical high state including high-energy 
INTEGRAL observations. On 5 January 2006, the trigger crite- 
rion was met when 3C 279 exceeded an R-band flux correspond- 
ing to R = 14.5 mag. INTEGRAL observations were carried 
out together with supplementing multifrequency observations 
from radio via IR and optical to X-ray energies. Preliminary re- 
sults from this campaign were reported by Coll mar et al.(2 007) 
Intensive ground-based monitoring in the IR and optical was 
assured by a WEBT (Whole Earth Blazar Telescope) cam- 



paign (see, e.g., Raiter i et al. (2006)[ |Villata et al. (2007)| and 
references therein). The results of the WEBT component of 
this campaign on 3C 279 in 2006 were presented in detail by 
Bottc her et al. (2007)[ In this paper we present details of the non- 
WEBT results of the radio, the X-ray, and soft y-ray bands (Sect. 
13, the interpretation of the radio data (Sect. [3]), as well as the fi- 
nal SED (Sect.HJl. In addition, we present a model interpretation 
of the SED, and discuss our results in the framework of leptonic 
emission scenarios (Sect.[5]l- 

Shortly after the core period of the multiwavelength cam- 
paign described here, the Major Atmospheric Gamma-Ray 
Imaging Cherenkov Telescope (MAGIC) reported the detection 
of very-high-energy (VHE: > 200 GeV) y-ray emission from 
3C 279 on 23 February 2006 ( |Albert et al. (2008)| l. This detec- 
tion was extremely surprising for two reasons: (1) FSRQs like 
3C 279 are known to have the peaks of their broad emission com- 
ponents in their SEDs at relatively low energies. In particular, the 
y-ray component of 3C 279 has been known from the EGRET 
era to peak at a few GeV, and the y-ray spectrum was generally 
not expected to extend out to VHE frequencies at a substantial 
flux level. (2) 3C 279 is the VHE y-ray source with the high- 
est confirmed redshift to date. At such large redshifts, yy ab- 
sorption of VHE photons by the extragalactic background light 
(e.g., |Dwek & Krennrich (2005)| l is expected to lead to severe 
attenuation of the spectrum in the VHE regime. This prompted 
Albert et al. (2008) to propose that the universe may be more 
transparent to VHE y-rays than previously thought (however, 
|Stecker & Scully (2 009) dispute this claim). It has been shown 
by |Bottcher et al. (2009)| that if this VHE detection of 3C 279 
can be confirmed, it would pose severe problems for any single- 
zone leptonic model interpretation and may potentially favor 
hadronic over leptonic models. 

However, the reported VHE y-ray flux, corrected for inter- 
galactic yy absorption, lies substantially above a straight ex- 
trapolation of non-simultaneous MeV - GeV spectra measured 
by EGRET as well as more recently by the Fermi Gamma-Ray 
Space Telescope ( Abdo et al. (2009)). Also, in spite of intensive 
efforts, no other VHE y-ray observatory has so far been able to 
confirm the detection of VHE y-ray emission from 3C 279. This 
indicates that the detection of 23 February 2006 must correspond 
to an extraordinary flaring activity and cannot be considered rep- 
resentative of the time of the simultaneous multiwavelength ob- 
servations reported here. Therefore, we will focus our modeling 
efforts on leptonic models, which have been shown to provide 
a successful description of all multiwavelength SEDs of 3C 279 
throughout the CGRO era ( |Hartman et al. (200 la")) . 

For the purpose of the model interpretation, we will 
use the widely accepted ACDM cosmology with Ho = 
70 km s _1 Mpc , Q m = 0.3, and £2a = 0.7. In this cosmology, 
the luminosity distance of 3C 279 is Di = 3.08 Gpc. Throughout 
the paper we refer to the energy spectral index a, defined through 
F v [Jy] oc v~ a , and the photon spectral index F p h = a + 1. 

2. Observations, data reduction, and results 

After 3C 279 met the trigger criterion for the targeted optical 
high-state (R < 14.5 mag) on 5 January 2006, we requested the 
pre-approved high-energy ToO observations with INTEGRAL 
and Chandra. To monitor the X-ray and UV-state of 3C 279 dur- 
ing the INTEGRAL observations, we proposed and were granted 
monitoring observations by Swift. In addition, ground-based 
monitoring observations in the framework of a WEBT campaign 
as well as sub-millimeter and radio monitoring of several tele- 
scopes were contributed to our campaign. In this section, we 



W. Collmar et al.: The Multifrequency Campaign on 3C 279 in January 2006 



3 



describe details of the collection and analysis of the previously 
unpublished data and briefly summarize the relevant aspects of 
campaign contributions which were previously published else- 
where ( [Bottcher et al. (2007)| |Bottcher et al. (2009)| l . 



analyzed as well (see Bottacini et al. (2007) ). We do not include 
them in this paper, because during the campaign superior optical 
measurements were performed by several ground-based obser- 
vatories (see Sect.ES) as well as the Swift UVOT (Sect. l2~3l 



2.1. INTEGRAL 

INTEGRAL is an ESA scientific mission dedicated to high- 
resolution spectroscopy with SPI ( |Vedrenne et al. (2003)| l and 
imaging with IBIS/ISGRI (angular resolution: 12' FWHM, 
point source location accuracy: - V - 3'; IBIS/ISGRI, see 
|Ubertini et al. (2003)| |Lebrun et al. (2003)} of celestial sources 
in hard X- and soft y-rays (15 keV to 10 MeV). In addition, 
INTEGRAL provides simultaneous monitoring at X-rays (3 - 
35 keV) with JEM-X ( |Lund et al. (2003)} and at optical wave - 
lengths (Johnson V-filter) with OMC ( [Mas-Hesse et al. (2003)] >. 

INTEGRAL observed the blazar in three pointings between 
13 January and 22 January 2006 for a total observation time 
of ~ 512 ks (Table [T). Due to annealing, no INTEGRAL 
SPI data are available. The INTEGRAL analyses were done 
with the INTEGRAL Offline Scientific Analysis (OSA) 
version 7.0. 





NGC 4593 


3C 279 


• 


□ 





Fig. 1. The ISGRI image shows a 6.4<x detection of 3C 279 in 
the 20-60 keV band. In addition, the Seyfert galaxy NGC 4593 
is even more clearly detected. The image shown covers a sky 
region of -11.8° x 5.6°. 



The IBSI/ISGRI data were screened for data failures or prob- 
lems, and in the end a total observation time of 505,043 seconds, 
available in 149 so-called science windows, was selected for fi- 
nal analysis. IBIS/ISGRI detected the blazar at energies between 
20 and 100 keV with a significance of ~1 .5cr during the total 
used exposure of 505 ks. The ISGRI image (Fig. Q]) of the 20- 
60 keV band shows a 6.4cr detection of 3C 279. Other Virgo re- 
gion sources are also detected in these observations, for example 
the prominent quasar 3C 273 and the Seyfert galaxy NGC 4593. 
The blazar was measured by ISGRI at a surprisingly low flux 
level of (2.53 ± 0.52) x 10" 4 ph cnr 2 s -1 . The spectral analy- 
sis with XSPEC between 20 and 100 keV, assuming a power- 
law shape, yields a weakly determined average spectral shape 
(Fig. [3]> of the photon index r ph = 2.0 + 0.5 (1(f). 

Subdividing the observations into three roughly equal time 
intervals (13 -15 January, 16-18 January, and 18-20 January), 
yields a low significance flux value for each interval. No signifi- 
cant variability is seen, but there appears to be a trend of a rising 
hard X-ray flux (20 - 60 keV) toward the end of INTEGRAL'S 
observational period (Fig. |4j». 

JEM-X did not detect 3C 279. Upper limits were derived for 
two energy bands (5-10 and 10-20 keV) based on the mosaic 
images from these observations. The OMC monitoring data were 




Fig. 2. The blazar is clearly visible in the Chandra image. A 
weaker source, probably the X-ray jet of 3C 279, is also vis- 
ible at an angular distance of ~4 arcsec. The purple (source) 
and the green (background) circles indicate the sky regions used 
for data analysis. The image shown covers a sky region of 
-34.0" x 24.5". 



2.2. Chandra 

After the INTEGRAL observations were approved, we also 
triggered our pre-approved Chandra 30 ks ToO observation on 
3C 279. It was scheduled on 17 January 2006, roughly centered 
in time on the one week INTEGRAL observation. Chandra ob- 
served the blazar for a total observation time of ~32 ks in order 
to determine contemporaneously and most accurately the X-ray 
state and spectrum of 3C 279. To avoid pile-up in the Chandra 
detectors of the assumed strong X-ray source, the LETG-ACIS- 
S mode was used. 

The blazar 3C 279 was significantly detected by Chandra. 
The Chandra X-ray image (Fig. I3 shows next to it a weaker 
X-ray source at a distance of ~4 arcsec. This could be emis- 
sion from its jet because it is along the direction of its VLA jet 



(e.g. Akujor et al. (1994) I. Chandra measured a well determined 
X-ray spectrum between 0.3 and 7 keV (Fig. [3j. Assuming the 
canonical power-law shape at X-ray energies, a spectral photon 
index of T p h = 1.55 + 0.02 was derived. The spectrum is consis- 
tent (reduced^ 2 : 1.14) with a hard power-law shape. However, 
the spectral analysis indicated a spectral bending from a harder 
to a softer spectrum toward higher energies. A fit with a broken 
powerlaw shape yielded a better fit (reduced^ 2 : 0.82): a harder 
shape (photon index: 1.46+0.03) below ~3 keV and softer one 
(photon index: 1.87+0.12) above. In particular, no soft excess 
was found, which would have indicated a contribution of the 
synchrotron component or the accretion disk to the soft X-ray 
emission. Chandra did not observe any variability of 3C 279 dur- 
ing the uninterrupted source-pointing of ~8 hours. 
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Fig. 3. X- to hard X-ray spectrum of 3C 279 as measured by 
Chandra (17 January 2006; exposure -32 ks), Swift (6 point- 
ings between 13 - 19 January 2006; exposure ~65 ks) and 
INTEGRAL/ISGRI (3 pointings between 13-22 January 2006; 
exposure ~505 ks) together with the Chandra (the Swift is com- 
pletely consistent with Chandra) and ISGRI best-fit power-law 
shapes (solid lines). The dashed line shows the extrapolation of 
the Chandra spectrum up to 100 keV, indicating the softer ISGRI 
spectrum above ~20 keV. The best-fitting photon indices includ- 
ing their errors (l<x) are given. 



2.3. Swift 

The y-ray burst explorer Swift was launched in 2004, 



and carries three co-aligned detector s (|Gehrels et al. (20 04)): 
the Burst Alert Telescope (BAT, Barthelm y et al. (2005) i, 
the X-Ray Telescope (XRT, |Burrows, et al. (2005)|> and th e 
Ultraviolet/Optical Telescope (UVOT, |Roming et al. (2005) I. 
Between 13 and 19 January 2006, Swift regularly observed the 
blazar, and — by our ToO request — contemporaneously with 
INTEGRAL (Table [TJ. The data analysis followed the standard 
procedure for Swift XRT data (e.g. |Grupe et al. (2010) 1 by using 
the HEASOFT 6.8 tools and XSPEC 12.5. l.n ( |Arnaud (T9 96)). 
All observations were performed in photon-counting mode 
fHill (2004)) l. Source photons were selected in the circular re- 
gion with a radius of 82.5 . Spectral data were re-binned to 50 
counts per bin (except for segment 003 where we used 20 counts 
per bin). 

The XRT measured a moderately variable 0.2 - 10 keV X- 
ray flux with a variability amplitude of AF/F < 25 %. After an 
initial decline throughout 13 January 2006, the Swift XRT flux 
seems to indicate a gradual rising trend during the remainder of 
the INTEGRAL observation period, in agreement with the trend 
seen by RXTE in the 2-10 keV range (see Fig. |4j). 

Parallel to XRT, UVOT monitored 3C 279's ultraviolet and 
optical emission. The UV magnitudes have been de-reddened 
using the extinction curves of Cardelli et a l. (1989)] with Ry = 
3.1 and A v = 0.095, taken from the NASA/IPAC Extragalactic 
Database (NECQ). The data indicate a continuous fading of the 
source in all six filter bands (Figs. [4] and 13), in agreement with 
the trend observed with ground-based optical observatories (see 
Sect.Ol. 
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Fig. 4. X-ray (top two panels) and optical/UV (lower panel) 
light curves of 3C 279 during the core campaign period, 13 - 
22 January 2006. Units for the X-ray light curves are: F2-10 
[erg cirT 2 s _1 ] and F20-6O [10~ 4 ph crrT 2 s~'] 



2.4. RXTE 

In the course of a long-term monitoring program led by A. 
Marscher of y-ray bright blazars using the Proportional Counter 
Array (PCA) on board the Rossi X-Ray Timing Explorer 
(RXTE), 3C 279 was observed with 2-3 pointings per 
week throughout (and beyond) the campaign period. Details 
of the observations and data analysis have been published in 
Chatterjee et al. (2008)| and the X-ray data have been discussed 
in the context of the WEBT campaign (see Sect. \2.5\ and the 
MAGIC detection in |Bottcher et al. (2009)[ The RXTE data are 
included in Fig. [4] and the multiwavelength light curves (radio 
through hard X-rays) throughout January and February 2006 are 
shown in Fig. [7] 

Before and during the core period of our multiwavelength 
campaign, the RXTE light curve showed 3C 279 to be in a low 
flux state, near its historical minimum, with a moderately soft 
spectrum with a ~ 1 in the 2-10 keV range. Toward the end 
of January 2006, i.e., after the core campaign period, the object 
went into a high-activity X-ray state with average X-ray flux a 
factor ~ 2 - 3 higher than during the previous low state. In this 
state, 3C 279 showed substantial, correlated flux and spectral 
variability on time scales of ~ 10 days. Notably, the time of the 
MAGIC detection occurred during the rising phase of a major 
X-ray outburst, peaking about one week thereafter at a flux level 
of ~ 5 times the previous quiescent flux. In the RXTE X-ray 
band 3C 279 is continuously variable with flux variations up to 
a factor of 10 (F2-10 between ~ 0.5 and 5 x 10~ n erg crrT 2 s _1 ). 
See e.g. |Marscher (2006)| for a longterm (1996 - 2005) RXTE 
X-ray light curve of 3C 279. 



1 http://nedwww.ipac.caltech.edu/ 



2.5. Optical monitoring 

The quasar 3C 279 is regularly monitored at optical, near- 
infrared, and radio frequencies by a number of observatories 
associated with the Whole Earth Blazar Telescope (WEBT). 
It was the monitoring provided by the WEBT optical partners 
that led to the discovery of the high state of 3C 279 on 5 
January 2006, which triggered our multiwavelength campaign. 
Early results of this campaign were published and discussed by 
Bottcher et al. (2007)] Here we add additional observations in 
the radio- to- sub-mm regime from radio observatories, which 
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Table 1. Observation log of the high-energy observations during our January 2006 campaign. For INTEGRAL the orbital revolution 
number is given. 



Obs. Periods 


T-Start 


T-Stop 


Exposure 


Satellite 


001 


2006-01-13 00:40 


2006-01-13 23:37 


15014 


Swift 


002 


2006-01-14 00:42 


2006-01-15 02:35 


10613 




003 


2006-01-16 04:07 


2006-01-16 23:34 


1581 




004 


2006-01-17 01:03 


2006-01-17 23:53 


19289 




005 


2006-01-18 01:14 


2006-01-18 24:00 


13897 




006 


2006-01-19 01:08 


2006-01-19 06:14 


4113 




001 


2006-01-17 08:59 


2006-01-17 17:51 


31920 


Chandra 


397 


2006-01-13 00:22 


2006-01-15 03:03 


174879 


INTEGRAL 


398/399 


2006-01-16 07:17 


2006-01-20 17:41 


322326 




400 


2006-01-22 05:13 


2006-01-22 09:25 


14392 





a Start and end times are given in UT 
b Exposure is given in s 



are not asscociated with the WEBT, and which complement the 
spectral and time coverage considerably. The combined radio, 
sub-millimeter, optical, and X-ray data are shown in Figs. |4]-|7] 

In the optical, the WEBT monitoring of December 2005 
- March 2006 indicated an overall elevated optical state, with 
variability amplitudes of ~ 0.5 mag on time scales of several 
days. The most remarkable light curve feature during that cam- 
paign was an unusually clean, quasi-exponential decay of the 
flux in all optical bands (BVRI) over a period of about two 
weeks (~ JD 2453744 - 2453757), including the core period 
illustrated in Figs. |4] and [7] shortly after exceeding the trig- 
ger threshold of R - 14.5 mag. This decay could be well 
approximated by an exponential F(t) - Foe~' /T with a time 
scale of t ~ 12.8 days, and continued throughout the time pe- 
riod of the high-energy (INTEGRAL, Chandra, Swift) observa- 
tions. [B5ttcheF&^incjpeT2009) suggest that this clean, quasi- 
exponential decay may be the signature of deceleration of a bal- 
listically moving plasmoid along the jet. 

The light curve recovered to a high level right after the period 
of the high-energy observations (around 22 January) and 3C 279 
remained optically bright until the end of March 2006. This in- 
cludes the time of the MAGIC detection (23 February 2006), 
during which the optical light curve did not show any remarkable 
variability features. While Bottche r et al. (20 07) found hints for 
spectral hysteresis and a hard lag within the optical bands in light 
curve segments later during the WEBT campaign, there was no 
evidence for optical spectral variability during the multiwave- 
length core campaign period (JD 2453748 - 2453758, i.e., 12 - 
22 January 2006). 

2.6. Radio to sub-millimeter variability 

The blazar 3C 279 was observed in the radio- to sub-mm bands 
by a number of different observatories. The data sets are com- 
piled from ongoing regular flux density monitoring programs 
and complemented by dedicated target-of-opportunity observa- 
tions performed during 2005 and 2007. The total of all flux den- 
sity measurements cover a frequency range from 1 .4 to 345 GHz, 
with the following contributions: The 100 m Effelsberg tele- 
scope of the Max-Planck-Institut fur Radioastronomie (Bonn, 
Germany) observed at 1.4, 2.7, 5, 8.4, 10.5, 15, 22, 32 and 
43 GHz, providing flux density measurements which are quasi- 



simultaneous in time (for details see: |Fuhrmann et al. (2 008) 



Angelakis et al. (2008) i. 



The IRAM 30 m telescope on Pico Veleta (Spain) con 
tributed with flux density measurements at 90 and 230 GHz 
These data resulted from the regular IRAM flux monitoring 
program (see Ungerechts et al. (1998) and Agudo et al. (2006) 
for a decription of the data reduction), and additional data ob- 
tained on the basis of target of opportunity observations pro- 
posed by us. These flux density measurements are comple- 
mented by additional measurements at 90 and 230 GHz with 
the 6x15 m IRAM interferometer on Plateau de Bure (France) 
and at 90 GHz with the Australia Telescope Compact Array 
(ATCA). To fill some gaps in the time coverage, we also included 
data from the VLA polarization monitoring, which is performed 
at 5 GHz, 8.4 GHz, 22 and 43 GHz ( [Taylor & Myers (2000)1 
We further included the 14.5 GHz data from the University of 
Michigan Radio Observatory (UMRAO) monitoring program 
and the 230 and 345 GHz data from the AGN monitoring of the 
Submillimeter Array (SMA) on Mauna Kea (Hawaii). The data 
from UMRAO and the SMA were already published in part by 
Chatterjee et al. (2008) and |Larionov et al. (2008)] where also 
the corresponding data reduction and calibration is described. 

In Fig.|5]we plot the flux density of 3C 279 in the radio bands 
versus time for the period 2004.4 - 2007.8. We have also added 



the optical R-band flux density from Chatterjee et al. (2008) In 
the lower panel of the figure we show the flux-density vari- 
ability of the VLBI core component at 43 GHz (same data as 



in |Chatterjee et al. (2008) but fitted with a circular Gaussian 
beam). Furthermore we add the times of zero separation from 
the VLBI core for the jet components C22 (2005.18 + 0.06), 
C23 (2006.415 + 0.15), and (2006.894 + 0.1 1) for C24 (see 
Chatterjee et al. (2008)||Larionov et al. (2008)| >. 



From Fig.|5jit is obvious that the variability amplitudes are 
more pronounced in the optical and mm-bands than in the cm- 
regime, which is consistent with a synchrotron self-absorbed ra- 
dio jet. Between 2005.9 and 2007.0 the optical to radio fluxes 
rose, indicative of an optical/radio flare, which peaked by the 
end of year 2007. During the rising phase of this flare, a first 
local flux density maximum at mm-wavelengths, respectively a 
plateau in the 14.5 GHz data was reached around 2006.4-2006.5, 
about half a year before the absolute maximum of the flux den- 
sities were reached (peak: 2006.9 - 2007.0). This local max- 
imum is clearly seen in the SMA data at 345 GHz (0.8 mm) 
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Fig. 5. Top: Radio and optical flux density plotted versus time. 
Different symbols and colors denote different observing frequen- 
cies: 14.5 GHz black circles, 90 GHz red squares, 230 GHz green 
diamonds, 345 GHz blue triangles, and optical R-band (ampli- 
tudes given in [mJy]) yellow triangles. Black squares mark VLBI 
component ejection as indicated in the panel below. Bottom: 
43 GHz VLBI core flux density versus time. Black squares with 
horizontal error bars denote the linearly back-extrapolated ejec- 
tion times of the VLBI components C22, C23 and C24 (in order 
of their appearance). 



and 230 GHz (1 mm) and also in the VLBI core flux at 43 GHz 
(7 mm). This indicates the superposition and blending of two 
flares separated in time by about 0.4-0.5 yr. It is tempting to asso- 
ciate the ejection of the VLBI components C23 and C24 with the 
observed double-peak structure between 2006 and 2007 in the 
total flux densities and in the VLBI core flux. After 2007.0, the 
optical flux density decayed very fast, followed by successively 
more shallow decays between sub-mm and cm- wavelengths. 

Between frequencies of 345 GHz and 90 GHz, the variabil- 
ity is fast and the light curves are partially undersampled in 
time. With the present time sampling, the fastest significant vari- 
ations appear on time scales of < 20 - 40 days. Faster varia- 




5.0 1 ' ' 

1 10 100 1000 
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Fig. 6. Quasi-simultaneous broad-band radio spectra of 3C 279. 
Different symbols and colors denote different observing epochs, 
as indicated by the legend box in the top left of the figure. 
Because of the irregular time sampling at the different observato- 
ries, the binning intervals differ from epoch to epoch. The spec- 
trum of November 2006 was observed near the time of the peak 
of the 2007.0 radio flare, see Fig. [5] The solid black arrow marks 
the spectral flaring component, which may be associated with 
the new VLBI jet components C23 in Larionov et al. (2008) 



tions cannot be excluded, but at this point remain speculative. At 
90 GHz a flux density increase from 12.8 in 2006.036 to 17.7 Jy 
in 2006.125 is observed in 34 days. This corresponds to a rate of 
0.14 Jy/day for the flux density scale. This in turn leads to an ap- 
parent brightness temperature of the variable component of Tb > 
8.4 x 10 12 K, which could be reduced to values below the inverse- 
Compton limit of 10 12 K, if the Doppler-factor is larger than 6 > 
4.4 (adopting the equations given in |Fuhrmann et al. (2~0 08)). 
With this Doppler-factor the corresponding linear dimension of 
the emission region would be 42/<5 = 9 .5 ^as, which corresponds 
to a linear scale of r < 7.9 x 10 17 /<5 = 1.8 x 10 17 cm. Using the 
above Doppler factor and a typical speed of the VLBI compo- 
nents of fi a pp 15 c, we obtain a minimum bulk Lorentz-factor 
of the jet of y m - m > 15 and an inclination of the inner jet with 
respect to the line of sight at an angle of i < 7°. 

2. 7. Quasi-simultaneous radio spectra 

From the available radio data presented in the previous sec- 
tion, it is possible to construct quasi-simultaneous broad-band 
radio spectra. Because of the irregular time sampling at the dif- 
ferent observatories, the binning of the data in time is differ- 
ent for each epoch and was chosen to obtain the best possi- 
ble frequency coverage. In Fig. [6] we show some examples of 
the radio spectra, which cover the frequency range from 1 .4 to 
350 GHz. The legend at the top left of the figure shows the ob- 
serving date and binning interval. The spectra are remarkably 
complex and show spectral evolution on time scales of weeks 
to months. The spectrum of December 2005 (circles) was ob- 
served just at the time before the 2006-2007 activity started and 
the 2007 flare began. The radio spectrum is flat up to about 
43 GHz, and steepens toward shorter (sub-mm) wavelengths 
(q , 43/225c/y.- - 0.25, a 2 25/345CH Z - 0.9. From January 2006 to 
May 2006, the spectrum becomes more and more inverted, with 
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Fig. 7. Multiwavelength light curves of 3C 279 during January and February 2006. The gray-shaded area indicates the core campaign 
period of the INTEGRAL, Chandra, and Swift observations; the vertical dashed line indicates the time of the MAGIC VHE y-ray 
detection. Units of the X-ray fluxes are: fVio f er S cm~ 2 s~'] and F20-6O [10~ 4 ph crrT 2 s -1 ]. The top panel shows the RXTE 2 - 
10 keV spectral index. The RXTE data are from Chatterjee et al. (2008) updated by A. Marscher (private comm.). 



a prominent spectral maximum near 100 GHz. In the context of 
the three-stage shock in jet model of |Marscher & G ear (1985) 
(see also Lindf ors et al. (2006)| l, this behavior can be interpreted 
as a spectral component, which is in its inverse-Compton stage 
(rising turnover flux S „, and decreasing turnover frequency y m ). 
Between 10 - 100 GHz, the observed spectral slopes vary con- 
siderably in all epochs. This suggests another variable spectral 
component with lower turnover frequency, which is less promi- 
nent than the one discussed before. In February 2006, this spec- 
tral component is best visible, causing a spectral break near 
43 GHz, and it may have propagated down to 15 GHz by May 
2006. Owing to the lack of multi-frequency VLBI data during 
this time, it is not possible to decompose the radio spectrum 
of the total flux into several spectra of distinct VLBI jet com- 
ponents. However, with regard to the published 43 GHz VLBI 



maps and the known jet kinematics (Chatterjee et al. (2008) 



|Larionov et al. (2008)[ ), it is very likely that the more promi 
nent spectral break near 100 GHz is associated with a new jet 
component, which must have emerged from the self-absorbed 
VLBI core region between mid-January and mid-February 2006. 
A likely candidate for an association of that spectral feature with 
a new jet component would be C23, which became visible at 
43 GHz after f = 2006.415 + 0.15 (see solid black arrow in 
Fig.|6]l. We note that an association of a secondary and less pro- 
nounced spectral break in the 20-30 GHz range with an evolving 



older jet component (C22 ejected at 2005.2) is a more specu- 
lative, but not unreasonable interpretation. 



3. Interpretation of the radio data 

The close correlation of the VLBI core flux with the total flux 
(see Fig. |5] and |Chatterjee et al. (2008)| |Larionov et al. (2008)| ) 
allows us to attribute most of the variability to be located within 
the synchrotron self-absorbed VLBI core, for which an upper 
limit to its size 6 is available {6 > 6Dil{ct var )l( \ + z), with lumi- 
nosity distance Di and variability time scale t var ): 6 < 0.03 mas 
or r < 5.7 x 10 17 cm. Applying the light travel time argument 
and following Jors tad et al. (2005)1 a lower limit to the rela- 
tivistic Doppler factor of 5 > 5.1 is obtained, which is close 
to the Doppler factors estimated in the previous sections. For 
synchrotron self-absorption the magnetic field can be calcu- 
lated via B ssa [G] « 2.3 x IOTAS' 2 ^6 /(I + z). This leads to 
B ssa ^15 mG for the VLBI core and inner jet, where we adopted 
v m = 90 GHz for the spectral turnover frequency (see previous 
section and Fig. |6). If we assume that most of the synchrotron 
energy is radiated near the turnover frequency v m , we can esti- 
mate the Lorentz factor of the electrons via the relation v,„ 



1 .2 x 10^ By 1 (B in [G], v m in [GHz]) to be y ~ 2200, which 
corresponds to an electron energy of E e = ym e c 2 1.1 GeV 
First order inverse-Compton scattering would boost this radia- 
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tion to Vi C K y 2 v m ~ 10 17 6 Hz, thus into the X-ray domain. This 
suggests that the X-ray emission of 3C 279 is dominated by 
synchrotron-self-Compton emission, while a different compo- 
nent (either Compton scattering of external photons or second- 
order Compton scattering of synchrotron photons) is needed to 
reproduce the gamma-ray emission near ~ 10 24 Hz, in agreement 
with our SED modeling (see Sect. |5). 



4. Spectral energy distribution 

The best multiwavelength coverage at IR/optical/UV frequen- 
cies was obtained on 15 January 2006, where observations in the 
near-IR, optical (UBVRI) and UV (UVW1, UVM2, UVW2) fre- 
quencies were available within ~ 1/2 hr of UT 05:00. We there- 
fore focus our attention on the simultaneous multiwavelength 
SED centered on that time. For the X-ray regime, we included 
the high-quality Chandra spectrum of 17 January 2006 in our 
SED, which is justified by the moderate X-ray variability dur- 
ing the core campaign. Given the substantial integration time re- 
quired to obtain spectral information from the INTEGRAL ob- 
servations, we use the time-averaged spectrum over the entire 
observation for inclusion in the SED. Justified by the long vari- 
ability time scales in the radio, data taken within a few days of 
15 January 2006 were included to extend the SED into the radio 
regime. The resulting SED, including radio, IR, optical (WEBT), 
UV (Swift UVOT) and X-ray (Chandra, INTEGRAL) and soft 
y-ray (INTEGRAL) is displayed in Fig. [8] and compared to the 
very-low-state SED from the multiwavelength campaign during 
INTEGRAL AO-1 (June 2003, |Collmar et al. (2004)) . 
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Fig. 8. Spectral energy distribution of 3C 279, centered on 15 
January 2006 (red downward triangles). All IR/optical/UV data 
were taken within +1/2 hr of UT 05:00. Also shown for com- 
parison is the SED from the INTEGRAL AO-1 campaign from 
June 2003 (Collma r et al. (2004)} . Despite the huge difference in 
optical flux, the X- and hard X-ray fluxes are surprisingly simi- 
lar. 



The most remarkable feature of the 15 January 2006 SED 
is the low X-ray flux, while the optical flux was clearly in an 
elevated state. The X-ray spectrum is virtually indistinguishable 



from the historical low-state SED of June 2003, while the optical 
flux is about a factor of 10 higher. The only slight difference in 
the X-ray spectral shape is the indication of a steepening toward 
lower frequencies, possibly indicating a low-energy cutoff of the 
X-ray spectrum (and hence, the underlying electron distribution, 
see Sect. [5}, for which there was no evidence in the June 2003 
SED. A comparison to a further historical SED is shown in Fig. [9] 

B 

Unfortunately, because the time of this campaign was before 
the launch of AGILE and Fermi, there was no simultaneous cov- 
erage of the critical MeV - GeV y-ray regime, where the high- 
energy component of 3C 279 is expected to peak. Therefore, 
constraints on the peak of the high-energy component can only 
be based on plausibility arguments in comparison with historical 
SEDs, an example of which is shown in Fig. [9] 



5. SED modeling 



As demonstrated in Har tman et al. (2001a)| a leptonic one-zone 
jet model, including both SSC and external-Compton scattering 
as high-energy radiation mechanisms, has been appropriate and 
very successful for modeling all simultaneous SEDs of 3C 279 
collected during the lifetime of EGRET. We therefore choose the 
same approach to interpret the low X-ray flux states presented in 
Fig. [8] including the SED from 15 January 2006, the center of 
the core period of our multiwavelength campaign described here. 
For this purpose, we use the equilibrium version of the leptonic 
SSC + EC model of |B6ttcher & Chiang (2002)| 

In this model, a relativistic electron population is assumed 
to be in an equilibrium state between on-going injection, escape, 
and radiative cooling in a spherical emission region of radius 
R, moving relativistically with bulk Lorentz factor F along the 
jet. Electrons are injected with a power-law distribution speci- 
fied by an injection spectral index q and low- and high-energy 
cutoffs j\ and j2, respectively. The radiative cooling is evalu- 
ated self-consistently with the radiative output via synchrotron, 
SSC, and external-Compton emission. Particles are assumed to 
escape on an energy-independent time scale f esc = r/R/c, where 
i] > 1 is a free parameter. Particle escape will become relevant 
for particles with energies below a critical break energy yi, de- 
fined through the condition that f esc = f coo i, the radiative cooling 
time scale. Escaping particles are leaving the emission region in 
random directions and will therefore primarily end up in a slower 
sheath that may be present around the fast spine producing the 
high-energy emission. The magnetic field B is a free parameter 
of the model, as are the parameters characterizing the external 
radiation field. For any given set of parameters, the code also 
evaluates the ratios of co-moving energy densities in the mag- 
netic field, u' B and the equilibrium electron distribution, u' e , the 
corresponding kinetic (L e ) and Poynting flux (Lb) powers, and 
their ratio, eg = u' B /u' e = LbIL 6 . For a more detailed description 
of this equilibrium model see |Acciari e t al. (2009)| 

The model fits only the near-infrared to gamma-ray emis- 
sion, which is believed to be produced at the sub-pc scales of 
the jet, where the source is still optically thick at radio wave- 
lengths. The optically thick radio emission is believed to be pro- 
duced on > pc scales. Above ~ 100 GHz however, the emission 
comes from sub-parsec scales of a region which is located at or 
very near the base of the VLBI jet (size of VLBI core at 100 
GHz < 10 18 cm). We find that the optical through X-ray SED 
of 15 January 2006 can be well represented with the following 



2 The interested reader may contact W. Collmar for 
an ASCII table of the figure. 
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Fig. 9. Fits to the two low-X-ray state SEDs included in Fig. [8j plus a historical low-X-ray state SED from the CGRO era in 
December 1992 / January 1993. The individual model radiation components for the 15 January 2006 SED are shown. See text for 
details. 



model parameters: L e = 1.2 x 10 45 ergs s l ,y\ = 1.5 x 10 3 , 
J2 = 10 5 , q = 3.7, B = 0.8 G, R = 3.5 x 10 16 cm, F = 15, 
and an external radiation field energy density w ext = 2.2 X 
10~ 4 ergs crrT 3 . In order to reduce the number of free param- 
eters, we choose the observing angle in a way that the Doppler 
factor D = (T[l -By cos O bs]) = T. The magnetic field spec- 
ified above corresponds to a value slightly below equipartition 
with eg = 0.56. The fit is illustrated by the red solid curve in 
Fig. [9] While the X-ray spectrum is always produced through 
SSC emission, the y-ray emission is strongly dominated by the 
EC emission. The model does, indeed, slightly over-predict the 
soft y-ray (ISGRI) flux. However, the ISGRI error bars are large, 
and our fit can still be considered marginally consistent with the 
ISGRI spectrum. 

Remarkably, the other two low-X-ray states can be modeled 
with essentially identical parameters, changing only the low- 



energy cutoff of the electron distribution. The fits to the P2 and 
June 2003 spectra shown in Fig. [9] have been generated choos- 
ing yi = 750 for P2 and j\ = 550 for June 2003. As already 
alluded to before, the higher value of the low-energy cutoff may 
be responsible for the break in the X-ray spectrum seen in the 
Chandra data of our January 2006 campaign. 

The change in y\ affects the peak frequencies of all radia- 
tion components as v sy <x y 2 , y EC oc y 2 , and vssc K j\, as long 
as Compton scattering occurs in the Thomson regime near the 
Compton peaks, which is the case for the parameters adopted 
here. In addition, increasing j\ also has the effect of increasing 
the synchrotron radiation energy density, hence leading to an 
increasing SSC contribution in the high-energy SED. In Fig. [9] 
we included the individual radiation components - synchrotron, 
SSC, EC (BLR), EC (disk), and the disk emission as dashed, 
dot-dashed, double-dot-dashed, double-dash-dotted and dotted 
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curves for the 15 January 2006 SED fit. This illustrates the dom- 
inant SSC contribution to the soft y-ray regime in this fit. In the 
states fitted with lower values of j\, the SSC component still 
dominates in the X-ray band, but plays virtually no role in the 
production of the y-ray emission. 

In the model calculations for 3C 279, all our cases are in 
the fast-cooling regime. In the fast-cooling regime, particles will 
cool below yi to a cutoff energy determined by a balance be- 
tween radiative cooling and escape. In an internal shock sce- 
nario, the value of j\ may be related to the relative Lorentz factor 
of two colliding shells, and the efficiency of transferring swept- 
up proton energy to accelerated electrons. Because we do not 
fully understand the details of this transfer process, this latter 
effect is commonly parameterized through electron acceleration 
parameters e e and both assumed to be on the order of ~ 10 _I 
(e.g., B6ttcher& Dermer (2010) ). Those parameters are defined 
so that e e is the fraction of swept-up proton energy transferred 
to relativistic electrons, and ( e the fraction of swept-up electrons 
that are accelerated to relativistic energies. Specifically, if the 
electron acceleration process results in an electron spectral in- 
dex of q > 2, the low-energy cutoff is given by 

m p q — 2 e e 

7i = r — (r shoc k - 1), (1) 

m e q — 1 £ e 

where r s h oc k is the Lorentz factor of the (internal) forward or re- 
verse shock resulting from the collision of two shells of relativis- 
tic ejecta, measured in the co-moving frame of the shocked ma- 
terial (Bottcher & Dermer (2010)| l. While r s h 0c k is a strong func- 
tion of the relative Lorentz factors of the colliding shells, e e and 
are expected to depend on the efficiency of the generation of 
turbulent magnetic fields mediating the energy transfer between 
protons and electrons behind the shock fronts. 

In |Bottcher et al. (2007)[ constraints on jet parameters could 
be derived from the estimated synchrotron peak flux in the 
15 January 2006 SED, and a hint for hard lags among the 
optical (BVR) bands. The synchrotron peak flux was trans- 
lated into a magnetic-field estimate according to Eq. 4 in 
Bottc her et al. (2007)| Inserting the values of our fit results, D = 
15 and eg = 0.56, that estimate becomes 

B>B es = 0.14(^ d )~ 1 G, (2) 

which perfectly agrees with the adopted value B = 0.8 G in our 
fit. 

Based on evidence for hard lags (R band leading the V and B 
band fluxes), Bottc heFet al. (2007)| developed even more param- 
eter constraints based on the assumption that the lags are caused 
by a slow acceleration process, which needs to be faster than the 
radiative cooling time scales. This led to the magnetic-field es- 
timate in Eq. 6 in Bottc heFet al. (2007)[ Using our fit values as 
well as the best-fit B - R time lag, Tbr = 3.75 d, we find 

B < B dcc = 0.20 (1 + ky 2/3 G, (3) 

where k is the ratio of Compton-to-synchrotron cooling rates 
(i.e., k = u' rdA lu' B if Compton scattering occurs in the Thomson 
regime). This upper limit is substantially lower than our adopted 
fit value. However, we point out that evidence for hard lags in 
the WEBT campaign data seems to have been most prominent 
in the second half of the campaign, where even hints for spectral 
hysteresis among the optical bands were found. No indication of 
spectral hysteresis was present during the time corresponding to 
the 15 January SED fitted in this paper. 



6. Summary and conclusions 

We have observed the prominent FSRQ 3C 279 over a pe- 
riod of about 10 days in an intensive multiwavelength cam- 
paign in January 2006. We obtained hard X-ray / soft y-ray cov- 
erage with INTEGRAL, additional X-ray monitoring observa- 
tions with Swift XRT, and RXTE, a high-quality 0.3 - 7 keV 
X-ray spectrum from Chandra, optical/UV monitoring with 
Swift/UVOT, and also long-term optical/near-IR/radio monitor- 
ing by the WEBT, and the Effelsberg -, UMRAO, IRAM, and 
ATCA radio telescopes, and the SMA. 

Our campaign caught the object initially in a high optical 
state, followed by a clean, quasi-exponential decline throughout 
the core campaign period. The X-ray monitoring observations 
indicate the opposite trend, namely a gradual flux increase. The 
optical and X-ray fluxes of 3C 279 clearly do not follow a corre- 
lated variability trend during our campaign. 

The campaign results allowed us to construct a well sam- 
pled SED from radio through soft y-rays. Comparing this SED 
to a historical SED from the EGRET era as well as INTEGRAL 
AO-1 (June 2003), we find that among EGRET P2 (December 
1992 - January 1993), INTEGRAL AO-1 (June 2003) and our 
January 2006 campaign, the optical flux level — but not the op- 
tical continuum spectral shape — had changed by almost one 
order of magnitude, while the X-ray flux and spectral shape ap- 
pear almost unchanged. Our January 2006 SED represents the 
highest optical flux state among the three SEDs compared here, 
and is the only SED in which we found evidence for a downward 
bending in the X-ray spectrum, possibly indicating a low-energy 
cut-off. 

We modeled the SED of January 2006 together with the 
other two low-X-ray flux SEDs mentioned above, with a one- 
zone leptonic model, including synchrotron, SSC, and external- 
Compton emission and self-consistently evaluating an equilib- 
rium electron distribution with the radiative output. This model 
was previously successfully applied to all simultaneous SEDs of 
3C 279 collected during the CGRO era. We find that all three 
low-X-ray SEDs can be well represented with this model. Most 
remarkably, the optical variability among these three SEDs can 
be reproduced by changing only the low-energy cutoff of the 
electron distribution, leaving all other parameters fixed. We cau- 
tion that this interpretation of the variability is obviously not 
unique, given the substantial number of parameters entering the 
model. However, the evidence for a low-energy break in the 
January 2006 X-ray spectrum lends additional support to the 
presence of a low-energy cutoff in the electron spectrum at ul- 
trarelativistic energies. This result may have important conse- 
quences for the evaluation of the total kinetic power carried by 
the jet in 3C 279. Given the steep spectral index of the elec- 
tron spectrum inferred from the steep optical spectrum, most 
of the kinetic energy in the leptonic component of the sub- 
pc scale jet is carried by the low-energy end of the electron 
distribution. Because we lack evidence (e.g., a bulk Compton 
component, Begelman & Sikora (1987)) for a substantial ther- 
mal plasma component in the jet, our results then suggest that 
the kinetic luminosity in relativistic electrons required for our 
model fit (L e = 1.2 x 10 45 erg s _1 ) provides a realistic estimate 
of the total kinetic luminosity in leptons in the relativistic jet. 
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